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We examine the possibility of recovering the supersymmetric (and soft supersymme- 
Oh, try breaking) Lagrangian parameters as direct analytical expressions of appropriate 
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physical masses, for the unconstrained (but CP and R-parity conserving) minimal 
supersymmetric standard model. We concentrate mainly on the algebraically non- 
trivial "inversion" for the ino parameters, and obtain, for given values of tan f3, 
simple analytical expressions for the fx, Mi and M2 parameters in terms of three 
arbitrary input physical masses, namely either two chargino and one neutralino 
masses, or alternatively one chargino and two neutralino masses. We illustrate and 
discuss in detail the possible occurrence of ambiguities in this reconstruction. The 
dependence of the resulting ino Lagrangian parameters upon physical masses is il- 
lustrated, and some simple generic behaviour uncovered in this way. We finally 
briefly sketch generalizing such an inversion to the full set of MSSM Lagrangian 
parameters. 
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1 Introduction 



In the Minimal Supersymmetric extension of the Standard Model (MSSM) JT], |2|], with- 
out additional theoretical assumptions, the necessary breakdown of supersymmetry leads 
to a large number of arbitrary parameters, consisting of all possible renormalizable soft 
supersymmetry breaking terms 0]. The latter arbitrariness simply reflects our present 
ignorance of the fundamental mechanism underlying the breakdown of supersymmetry. 
Although several theoretically appealing scenarios have been explored, with additional 
assumptions considerably reducing the arbitrariness of the soft SUSY-breaking sector 
(typically like in hidden sector supergravity models |§, or gauge-mediated susy break- 
ing models ||), alternative strategies may be useful to exploit the information from the 
present and future collider data. If SUSY is realized at low energy in some way, and some 
super-partners discovered, the next immediate task would be to reconstruct from the data 
the structure of the SUSY and soft-SUSY breaking Lagrangian [that is, to determine as 
precisely as possible how the experimentally measured parameters would translate as pos- 
sible values (or bounds) on the MSSM Lagrangian parameters.] 

The expressions of physical parameters (mass eigenvalues, mixing angles and physical 
couplings) as functions of the most general, unconstrained MSSM Lagrangian, are well 
known BJ. Now, inverting those relations (i.e. expressing the Lagrangian parameters 
in terms of physical parameters) is non trivial in the unconstrained MSSM, due to the 
large number of parameters and possible mixing among them. This is especially true in 
the ino sector which involves a relatively complicated structure of the mixing terms, with 
e.g. a 4 x 4 mass matrix to "de-diagonalize" in the neutralino sector. To our knowledge, 
no systematic analytical inversion has been explored up to now. 

The aim of the present paper is twofold: 

• Investigate as much as possible analytically the reconstruction of the Lagrangian 
parameters from the physical masses, especially in the ino sector, identifying clearly 
the procedure and the related ambiguities when they occur; 

• construct a numerical code based on algorithms using as much as possible the above 
mentioned results, and which can be readily used for a full numerical study of the 
Lagrangian parameters as a function of the physical masses. 

For defmiteness we will restrict ourselves throughout the paper to the (CP and Im- 
parity conserving) unconstrained MSSM. Also to simplify the presentation we will con- 
sider two illustrative "scenarios" where either two chargino and one neutralino masses, 
or one chargino and two neutralino masses are input. The outcome is a well-defined 
algorithm providing, for given tan/3, the values of the /i, Mi and M2 parameters in terms 
of three arbitrary input masses. Furthermore, one of our main results is the derivation, 
in terms of simple analytic expressions, of the full neutralino mass spectrum and one of 
the Lagrangian parameters (Mi), when the three other parameters (M2, /x, tan f3) and one 
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neutralino mass are given. This together with the analytical expressions of the chargino 
sector constitute the building blocks of our algorithms. 

It should be stressed here that we are not aiming, in this paper, at a specific phe- 
nomenological reconstruction of the MSSM parameters from experimental observables. 
Clearly an algebraic approach cannot replace more systematic phenomenological studies 
of the reconstruction of the basic MSSM parameters. Various experimental strategies to 
extract the MSSM physical spectrum from present and future collider data have been 
considered, together with a detailed assessment on mass measurements, e.g., in ref. [|7|, H]. 
On the other hand, a complete scanning of the unconstrained MSSM parameters would 
be in practice rather tedious and extremely time consuming. We believe that our study 
would be useful and rather flexible for such reconstructions: besides the fact of provid- 
ing a fast numerical inversion routine, we systematically encode the possible ambiguities 
that can arise through the reconstruction from some of the physical masses. In the most 
optimistic case (that is, if knowing a sufficient number of physical input ino masses and 
with sufficient accuracy), this inversion allows for a precise reconstruction of the uncon- 
strained MSSM ino sector Lagrangian. If a more limited experimental information of the 
ino spectrum were available, such as only mass bounds or mass differences typically, one 
could still derive corresponding constraints on the MSSM Lagrangian parameters as well. 



Our input/output parameters are deliberately chosen so as to render the analytical 
inversion the most simple and transparent. As will be explained in details later on, our 
algorithms will be either fully analytical, in which case at most a readily solvable cubic 
equation is involved (in contrast with the generically quartic one in the ordinary diago- 
nalization of the neutralino sector or needs some numerical iteration, in which case 
it relies essentially on a system of two algebraic equations, that are only quadratic (in 
contrast with a "brute force" inversion, which would involve a highly non-linear set of 
equations). It thus demonstrates the feasibility and relative simplicity of such a system- 
atic inversion and more importantly we also exhibit the general trend and sensitivity of 
the Lagrangian parameters to the physical masses. Note, however, that we do not assume 
at this level any knowledge of the couplings of inos to higgses and gauge bosons. These 
would of course lead to extra information |T(| which together with our results could allow 
an even more complete reconstruction and cross-checks for the ino sector. 



The rest of the paper is organized as follows. In section 2, we briefly recall the main 
content of the MSSM Lagrangian, just to fix our definitions and conventions. In section 3, 
we present our main results for the analytic inversion of the ino parameters, and discuss 
uniqueness conditions. Some essentially technical material is given in two appendices. 
In section 4 we illustrate the dependence of the ino sector Lagrangian parameters, as 
functions of the relevant physical input masses. We also illustrate some resulting values 
of the ino Lagrangian parameters at a grand unification (GUT) scale. In section 5 we 
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briefly outline, for completeness, similar inversion relations in the scalar fermion and Higgs 
sector of the MSSM. Finally section 6 presents the conclusions and an outlook. 



2 MSSM Lagrangian and spectrum: reminder 

The supersymmetric part of the unconstrained MSSM Lagrangian consists of the gauged 
kinetic terms for the £77(3) x SU(2) x £7(1) vector multiplets, involving the three gauge 
couplings g 3 , g 2 , and gi, respectively, the gauge invariant kinetic term of matter fields, 
and the superpotential: 

W = uYtQHu - dY b QH d - eY T LH d + fiH u H d . (2. 1) 

In eq. (|2.1[ ), the hats denote superfields, Yi are the Yukawa couplings, and H u H d = 
tijH^H^ (with E\2 — 1 ) fixes our sign convention for fi. In (|2.1|) and in the following we 
will in fact neglect all flavor non-diagonal terms. We also suppress any flavor, color, etc 
indices. Moreover, we will restrict here to the case where all parameters are assumed to 
be real. 



The soft SUSY-breaking terms involve 
-the trilinear coupling terms: 

-{uA t Y t H u Q + dA b Y b H d Q + eA T Y T H d L + h.c.) ; (2.2) 

involving the Ai parameters (which have a mass dimension in this convention); 
-a contribution to the sfermion mass terms: 

Csfermion = ~M 2 Q (t* L t L + b* L b L ) - M? R t* R t R - M^ R b* R b R (2.3) 

-Ml{flf L + vlv L )-M 2 TR f R f R 

and the gaugino mass terms: 

Mi - M > - - M < - - 
C„ no = ~^BB - -fW'Wi - -fG a G a ■ (2.4) 

which fixes our sign conventions. There are also supersymmetric contributions to the ino 
masses, of which we write here for illustration only the ones contributing to neutralinos 
(after electroweak symmetry breaking): 



^neutralino = m z c w sin pW 3 H u - m z c w cos PW 3 H d 

+m z s w cos [3BH d — m z s w sin [3BH U + jj,H u H d 

(2.5) 

with s w = smO Wl c w = cos 0w, and tan/5 = v u /v d the ratio of the two Higgs vacuum 
expectation values, v Uyd = (H u ^ d ). It is now easy to see from Eqs.(^l, |2~3| ) that only the 
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relative phases of say (M 2 ,Mi) and (M 2 ,/i) are physically relevant. Indeed any phase 
change of M 2 in Eq. (|2.4| ) can always be absorbed by a phase change of the W field. The 
latter however fixes uniquely the phase change of H4, H u and B in Eq.( |2.5| ) in such a way 
that the phases of the combinations M 2 /Mi and [iM\ remain unchanged. The discussion 
is of course reminiscent of the one carried out for instance in 0] . In the present study, as 
already mentioned, we restrict ourselves to real-valued parameters, but do not necessarily 
assume universality of gaugino masses. We can thus choose, without loss of generality, 
M 2 to be always positive and Mi and fi to have arbitrary signs. 

Finally, the Higgs potential is built from soft SUSY breaking and F-term contributions 
to the Higgs scalar mass terms plus quadrilinear D-terms, and reads, before SU(2) x U(l) 
breaking: 



Vm 99 s = {m 2 Hd + ^)\H d \ 2 + {m 2 Hu + fi 2 )\H u \ 2 + {BfiH u H d + h.c.) (2.6) 

rl-h 
eff 



id 

1 9l + 92 /\ TJ |2 |tt|2\2 1 92 I tt« TJ |2 1 xrl-loop 

\ g K\ H d\ -Wu\) +—\±i d H Uji \ +V eff 



which also fixes a sign convention for the B parameter, and where K// °° P i s the one-loop 



contribution to the effective potential |TT|, 112, [T3|, The corresponding Higgs masses 



and mixing angles (after electroweak symmetry breaking) are given in Appendix 1. 



3 Inverting the ino MSSM spectrum 

Let us first sketch our general procedure to reconstruct the ino sector parameters. As 
already mentioned in the introduction, we have to fix a specific choice of input /output 
parameters, that we do as follows. We first assume that tan/3 is an input parameter at 
this stage, i.e. that it has been extracted from elsewhere prior to ino reconstruction, or 
simply fixed arbitrarily. (Obviously, once a reconstruction algorithm is defined for fixed 
tan/3, one may always study the sensitivity to this parameter). Then, we consider two 
basic scenarios, Si and S2, and the corresponding algorithms. 

Si 1 , the input are two chargino masses, M x +, M x +, and one (but any) of the neutralino 
masses, ±M^- 

S2: the input are a single chargino mass, M x +, and two neutralino masses, say ±M^ 2 , 
±M N3 . 

It is important to note here that we will adopt the formulation in which the eigenvalues of 
the chargino mass matrix are by construction always positive, while those of the neutralino 
mass matrix can have either signs. This means that in the inversion process one should 
consider both ±M^ i as possible inputs [a feature consistently taken into account in the 
formulae we derive.] 
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For each of the above scenarios the aim will be to identify a corresponding algorithm 
which allows the determination of the output values of fj,, M 2 , Mi, as well as the remaining 
three neutralino masses (resp. two neutralinos and one chargino masses) in the Si (resp. 
S 2 ) case. 

First, as we will show in detail below and in Appendix [FJ, scenario S% allows a fully 
analytical algorithm to handle all the features of the inversion procedure. Algorithm Si 
will therefore be our starting point for the whole study []. We stress, however, that it does 
not necessarily imply a strong prejudice on the "chronology" of discovery of the inos since, 
as we shall explain, it may also be used as a basic building block for a probably more 
phenomenologically relevant situation, in scenario 5*2 below. [Furthermore, we will see 
that Si can be naturally separated into two independent steps corresponding to chargino 
and neutralino sectors. It can thus also be readily used in a more general context than 
the one of the present paper]. 

To scenario S 2 we associate the following algorithm: first assuming temporarily that 
e.g. fi is an input parameter, together with two neutralino masses, say Mjv a and M^ 3 , 
then a simple quadratic system gives M 2 , Mi as function of the latter input, as well as 
all other ino masses. Now, the key observation is that it is relatively simple to merge 
the solutions of the latter system with the previous algorithm Si, to obtain consistently 
\i, M 2 , Mi from a single chargino mass, M x +, and two neutralino masses, say M^ 2 , M^ 3 
input. More precisely, choosing an adequate "initial guess" value for e.g. M 2 one can 
simply use Si (actually only a part of it) to determine /1 in terms of (M +, M 2 ), followed 
by the above mentioned Mi, M 2 solutions, iterating with respect to M 2 until a convergent 
set of values is obtained Q. While Si was fully analytical, the price to pay is that, accord- 
ingly, 5*2 has to be partly numerical through the iterative procedure. In most cases, this 
combined algorithm turns out to converge rapidly, after 2 to 3 iterations for an accuracy 
that is sufficient for all practical purposes. 

This particular decomposition, with this choice of input/output, is the one giving 
the most algebraically tractable inversion in the ino sector. In fact, while Si alone is 
probably not very relevant physically (since for a rather generic choice of the Lagrangian 
parameters, it looks more likely to have two neutralinos and one chargino in the lightest 
part of the ino spectrum |7|, Q, jl4l), precisely this situation is tractable from the combined 
algorithm of scenario S 2 , as explained above. However, it is rather instructive to study 
algorithm Si separately in some detail, as it exhibits important properties of the inversion 
in a relatively simple manner. [Of course it is also possible to use any part of our particular 

Except when ambiguous, we call from now on the scenario and its associated algorithm with the 
same name. 

2 There are always (complex) solutions since, as easily checked, the explicit non-linear system equivalent 
to S2 is a (high degree) polynomial. Non convergent domains thus simply correspond to the impossibility 
to match our basic assumption of real [A, Mi, M2. 
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procedures to simply scan over some values, if not known, of some of the ino masses, as 
will be illustrated in the plots in section 4]. 



3.1 Basic algorithm Si'. M x +, M x +, M N2 input 

Assuming tan j3 and the two chargino masses given, one easily obtains from the expressions 
of the chargino mass eigenvalues (see eq ( |A.2|) in Appendix |A]): 



/i 2 = i(M 2 ++M 2 +-2m 2 w ±[(M 2 ^ (3.1) 



and 



M 2 = [M 2 + + M\ t - 2m 2 w - /i 2 ] 1 ^ (3.2) 



In eqs. ( p.l|) and (|3.2|) the ± ambiguity in front of the square root results from the 
invariance of the physical masses (see eq. ( |A.2| ) of appendix A) under the substitution 



<-> M 2 (we stress again our sign convention, M 2 > in (|3.2j )). In other words, from 
the two chargino masses input only, there is clearly no way of assessing the amount of 
gaugino or higgsinos components of each of them. To achieve this, one would require 
the knowledge of mixing angles (i.e. couplings) |10[ or alternatively of two neutralino 
masses, as will be discussed below in scenario S2 ■ In what follows, we thus arbitrarily 
choose to illustrate only the case |/i| < M2 (corresponding to the minus sign choice in 
front of the square root of ( |3.1| )), i.e. the case where the lightest chargino has a stronger 
higgsino like component. We thus refer to this situation as "higgsino-like" . The resulting 
output of algorithm Si corresponding to the opposite "gaugino-like" situation, > M 2 , 
is therefore trivially obtained by interchanging the values of and M 2 (together with 
the correct sign of /i assignment, see e.g. captions of fig. Q). In addition, since M 2 > 0, 
Mi will have arbitrary sign (see previous section 2). [Had we chosen to let the sign of M 2 
rather than Mi undetermined at this stage, our algorithm would have been slightly more 
complicated] . 

Concerning now the sign ambiguity on /i, it is in fact fixed from another relation from 
eqs. ( |A.2|) (implicitly used in eq. (|3.1|)): 



M 2 ^ = m 2 w sin 2(3 ± M x +M x + (3.3) 

for each respective choice of the ± ambiguity. The latter ±, also appearing inside the 
square root in (|3.1| ), corresponds to a true ambiguity, i.e. whenever the expression under 
the square root of eq. ( |3.1|) is positive for both sign choice, there are two distinct solu- 
tions for (fi, M 2 ). Obviously, the occurrence of this discrete twofold ambiguity depends 
crucially on the values of M x + , M x + and tan (3, and deserves a more careful examination 
to which we now turn. 

We illustrate in table 1 and fig.[I] the different domains where one obtains either no 
solutions (for /1 real), a unique solution, or the discrete twofold ambiguity, depending 
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on the input values M y 
difference A v + = |M V + - M v , 



M + and tan/3. What is relevant here are the chargino mass 

^2 



and sum S x+ = M + + M Y +, as well as the quantity 



€ 1 2M Y+ M Y+ + € 2 (M 

X-l A2 



+ MK - 2Ml 

X-2 



W) 



2M^ 



(3.4) 



where e\,€2 = ±1. For instance, the last two columns of table 1 summarize respectively 
the occurrence of the twofold ambiguity or the existence of (real) solutions at all. In some 
zones this is fully controlled by the values of tan (5 which we expressed in the form of 
necessary and sufficient constraints on sin 2/3. This is systematically taken into account 
in our algorithm, which gives in the relevant case the two possible solutions for /i. One 
then needs an extra information to eventually eliminate one of the solutions. Note finally 
that in fig.|l] we assumed for simplicity that M x + > M x + which is obviously just a matter 
of convention. [In any case, the opposite situation corresponds simply to a symmetry 
around the bisectrix line.l 



A X+ /V2M W 


A X +/2M W 


S X+ /2M W 


X+<0;X-<1 
in all cases 


twofold 
ambiguity 


constraint 

on (3 


< 1 


< 1 


< 1 


0<X^<1 


No 


X-<sin2/3<Xr 


zone (I) 






X+ <0 




< 1 


< 1 


> 1 


o<x-<i<xr 


No 


sin 2(5 > X- 


zone (II) 






X±<0 




> 1 


< 1 


< 1 


exp. excluded 






> 1 


< 1 


> 1 


X~<0; 1<XZ 


Yes when 


No 


zone (III) 






Q<X+<1 


sin 2(5 < Xi 




> 1 


> 1 


> 1 


X-<0;l<Xl 


Yes 


No 


zone (IV) 






Xi > 1 







Table 1: Inequalities characterizing the different domains of solutions /i, M 2 for arbitrary 
M v+ , M v+ input. A y+ = |M v+ - M+l; S y+ = M v+ + M v+ . Xi is defined in the text. 

Xi 1 X2 * 1 X2 Xi " AT Xi X2 ^ 



3.2 Neutralino mass inversion 

Let us now turn to the de-diagonalization of the neutralino sector. The question we want 
to answer analytically here, is how to determine M\ and three neutralino masses, when 
M2, fi, tan/? and a fourth neutralino mass are given. We will only outline the procedure 
leaving some technical details to the appendix. As we shall see, the discussion in this 
section is quite general and serves indifferently as a basis for the more specific algorithm 
Si or S 2 . 
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M 2 M W 



M x 2 - M xr V2 Mw , 



M X2= 2 y M w 



| experimentally fully excluded region 



HI excluded if |i and M2 are both real valued 



M. 



Figure 1: Domains of solutions for //, M 2 as functions of M x +, M x + and tan/3, where we 
assumed without loss of generality M x + > M x +. The detailed definitions and properties 
of the different domains (I)-(IV) are explained in Table 1. 



Since we restrict ourselves to the case where Mi , M 2 and \i are all real-valued, the 
neutralino mass matrix 



M 



( Mi — m z sw cos ft mzSw sin ft \ 

M 2 m z cw cos ft —m z cw sin /3 

— m z sw cos /3 m z cvK cos ft — /i 

V m^s^y sin ft —m z cw sin /3 — /i / 



(3.5) 



is symmetric and can be diagonalized through a similarity transformation, i.e. 



PMP- 1 = M diagonal (3.6) 

Analytical results for this diagonalization are well-known Q. However, a straightforward 
inversion of the expressions of the mass eigenvalues using these results is far from obvious. 
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A more manageable alternative is to start from the four quantities 



TrM (3.7) 

(TrM) 2 Tr(M 2 ) , . 

(TrM) 3 TrM Tr(M 2 ) Tr(M 3 ) 

- 6 2 + (3 ' 9) 

DetM (3.10) 



which are invariant under any similarity transformation. These quantities allow one 
to relate the four mass eigenvalues to the initial parameters in ( |3.5|) as follows: 



M Nl + M Na + M Na + M Na = M 1 + M 2 (3.11) 
M Nl M N2 + M N2 M N3 + M N . A M Ni + M Ni M Nl = - // - M% (3.12) 



M Nl M N2 M Na + \I x ,.\l x ,.\I Xl + M N3 M Ni M Nl + M Ni M Nl M N2 

= /.Mf sin 2/3 - (n 2 + c 2 w M 2 z )M 1 - (/i 2 + s 2 w M 2 z )M 2 (3.13) 

M Ni M N2 M Ns M Na = iiM 2 z {c 2 w M x + s 2 w M 2 ) sin 2(3 - ii 2 M x M 2 (3.14) 



Here M^. denote the eigenvalues of the mass matrix Eq. (|3.5|) , to be distinguished from 
the physical neutralino masses Mat. given by = e^M^^ where = ± is undetermined 
at this level. 

We emphasize that Eqs. (|3.11| - |3~T^) give the complete information on the relationship 



between the four neutralino mass eigenvalues and the original parameters, //, Mi, M 2 and 
tan j3. Furthermore these equations are a good starting point for our purpose, as they do 
not favour a priori any particular set of parameters. Thus, the system may be solved in 
many different ways, depending on the choice of input/output one is interested in. 

According to the algorithm Si, described in the previous section, we can assume that 
jj,, M 2 (and tan/5) are determined at this stage, and extract Mi and the three physical 
masses M Nl , M N3 , M Ni as functions of one mass, say M N2 [] 

Accordingly, in this Si scenario, one finds for Mi 



3 Actually what we call M^ 2 plays the role of any neutralino mass to be given as input. There will be a 
relabeling of neutralino states depending on the values of other parameters, and one can check afterwards 
whether the input mass was the lightest, heaviest or intermediate one. 
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Mi 



Pii + P 2l (v 2 + M 2 Z + M 2 S 2l 



si) 



+ ^M 2 z M 2 s 2 w sin 2/3 



P2i(S- 



2i 



Mo 



+ H(c 2 w Mf sin 2(3 



HM 2 



where 



2 1 



M N2 + M 



Ni 



(3.15) 



P 2i = M N2 M Ni 

and Mtv, is any of the three remaining neutralino mass parameters. 
The latter take the following form, 



M 



1 



3ai 



(a 2 - 2 v /Z Acos(Ar^[ J B])) 



(3.16) 



M: 



7T\ 



3ai 



(a 2 + 2 v /z Acos(Ar^[5] - ^)) 



(3.17) 



M; 



JV 4 



3ai 



ja>2 



2V- 



71 , 



■Acos(Arg[B] + -)) 



(3.18) 



where A, B, a\ and a 2 are given in Appendix |B|, Eqs.( |B.7| , BTTI , |B.12| ), and Arg[B] 
denotes the phase of the generically complex valued quantity B. Note that M 1; M Nl , M Na 
and Mtv 4 are all real- valued since we chose /x and M 2 to be real. Note also that Eq.( |3.15| ) 
should yield the same value when M^ is substituted by any of the three neutralino masses 
given by Eqs. (gig - gig) . This allows a non-trivial consistency check. For a more detailed 
discussion of the derivation of Eqs.( 3.15 - 3.18 ) and related material the reader is referred 
to Appendix |B[ 



3.3 Combined algorithm scenario 5^: M x +, Mn 2 , Mjv 3 input 

Alternatively, when considering now scenario S 2 as defined previously (i.e. the inputs are 
fi, Mjv 2 , Mjy 3 ), an even simpler solution occurs for M 2 and Mi from eqs ( p.ll| )-( p.l4| ). In 
fact, the necessary and sufficient conditions for the existence of solutions to eq. ( p.llQ - 
( |3.14|) give the consistency conditions eqs ( B.4|) , ( |B.5|) (see Appendix Q), from which one 
immediately obtains, for given /j, simple quadratic equations for M 2 , Mi, that we omit 
to write here. The price to pay, as already mentioned previously, is that without further 
model-dependent assumption, this combined algorithm for scenario S 2 has to rely on a 
numerical (iterative) consistency check, and some non-trivial procedure has to be per- 
formed in order not to miss all possible solutions The upshot is that up to at most 
four distinct solutions for (//, M 1; M 2 ) can occur for given M x +, M N2 , M Ns input, once 



Indeed, solving e.g. eq. fl3.1| ) for /i as a function of M x +, M2, and bravely injecting the solution 
in eqs. (B.4)-(B.5) gives a highly non-linear equation for, e.g., M2. Most of the solutions are in fact 
spurious, and our iteration algorithm, based on two equations which are only quadratic, copes with these 
redundancies in a simpler way. 
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all consistency constraints (including our sign convention choice M 2 > 0) are taken into 
account. Unfortunately, in contrast with the twofold ambiguities of algorithm Si alone, it 
is not easy to work out analytically specific domains of the input parameter space corre- 
sponding to a definite number of distinct solutions, and one has to rely on the numerics. 

Let us summarize this section: 
-in scenario Si, we start from M X +,M X +,M^ 2 to determine first /i and M 2 (up to the 
possible twofold ambiguity), and then M\ from eq. ( |3.15| ), given tan/3. The solution for 



Mi, ( p. 15 ), is indeed unique for fixed M 2 , \i and one neutralino mass values. This part 



of the algorithm may thus be used in a more general context (than Si), where \i, M 2 
and tan/? would have been extracted in a way or the other. Within scenario Si alone, 
there are in principle two possible M\ values for any two chargino and one neutralino 
mass input, since the chargino mass input does not distinguish the gaugino-like from the 
higgsino-like situation, and ( |3.15|) is not symmetrical under the interchange |//| «-> M 2 . 



-Alternatively, in scenario S" 2 one can obtain /i, M 2 , M\ consistently by iteration from 
M x +, Mat 2 , and M^ z . In practice, this combined algorithm converges very rapidly, but 
a relatively large number of distinct solutions (up to four) can a priori occur in some 
domains of the relevant input parameter space (/z, Mi, M 2 ). [Note that these ambiguities 
are different from the previous Higgsino-like <-> gaugino-like ambiguity of scenario Si, the 
latter being precisely removed from the knowledge of a second neutralino mass]. 



4 Numerical illustrations of the ino inversion 

We shall illustrate here with some representative plots the results of the inversion in the 
ino sector, according to the algorithms explained in section 3, together with a few remarks 
and comments. Since our choice of input masses is rather arbitrary and may not directly 
correspond to the most interesting experimental situation, our comments are accordingly 
essentially qualitative. We will see, nevertheless, that a number of general and interesting 
properties of the inversion can directly be seen here, irrespective of the precise values of 
the other parameters that have to be fixed, like tan (3 typically. 



4.1 Scenario Si', two charginos plus one neutralino input 

We first illustrate the basic algorithm Si, namely the reconstruction of /i, M 2 (and Mi) 
from M x +, M x + (and tan/? and M N , 2 ) in fig. 0. As already mentioned, we illustrate only 
the higgsino-like situation \ < M 2 since the gaugino-like illustration is trivially obtained 
from it (see captions of fig. §). 

Since it is phenomenologically more likely [0] that the first inos discovered (if any) will 
be two neutralinos and only one chargino (the second chargino being of heavier mass), we 
fixed in scenario Si only one chargino mass, say M x +, while varying the other one, M x +, 
in a large range to illustrate as much as possible the dependence on the physical input 
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(see figure captions). The shape of the various plots in fig. 2 is quite generic. First, it 
exhibits three distinct zones as regards the existence, uniqueness, or possible ambiguities 
on /x, M 2 , Mi (see also the discussion and table 1 in section 3): 

i) the grey shaded region, where there are no solutions for real /x, corresponds to zone 
(I) of table 1, with the corresponding constraint on sin 2/3 not fulfilled (]. If one 
takes a smaller or larger M x +, this region around M x + will be simply displaced 
accordingly. 

ii) in the left and right boarder zones are the regions of twofold ambiguities on /x, M 2 , 
as indicated. Note therefore that one of the two solutions has a discontinuity at 
the boarder between the single and twofold solutions region. Without additional 
information (or particular model assumption) one cannot a priori reject any of the 
two solutions. 

iii) Finally the two bands in between zones i) and ii) correspond to the region where 
eqs. ( |3.1|, [3l?D give a unique solution for /x and M 2 . It is interesting to note that 



those bands are narrower when tan (3 is increasing (in fig. 2 tan/3 = 2), irrespective 
of the M + values, becoming e.g. only a few GeV wide for tan (3 > 35. 

Moreover, as a general behaviour, the values of /x and M 2 are rather insensitive to tan /3 
(apart obviously from the discontinuous change occurring for one of the solution at the 
boarder between zones ii) and iii)). One can also remark from fig. 2 the relatively 
simple shape of /x and M 2 as function of M x +, with an almost constant or linear depen- 
dence, apart from some narrow regions. This is easily understood, since from eqs. ( |3 . 1| ) , 
one obtains (p,M 2 ) ^ (M +, M +), (resp. (M +,M +) ) for M + « M + (resp. 



x 2 ' xi n ^ 1 v xi ' x 2 ' ' x 2 xi 

M xt » M x +). 

In fig. |2] we also plot Mi for the corresponding values of /x and M 2 , and for fixed 
Mjv 2 . The twofold valuedness of Mi there is entirely due to the twofold ambiguity of /x 
and M 2 . In contrast, we do not illustrate here the other determination of Mi resulting 
from the \fi\ «-> M 2 interchange: although the resulting Mi may be in general quite 
different, this interchange will have very little effects as a function of M x +. First of all, 
the trivial, almost constant behaviour of Mi in this plot (except for small M x + < 100 
GeV) is relatively simple to understand from eq. ( |3.15| ). Indeed, M 2 is large for large M + 
(since we assume |/x| < M 2 ). In this limit and for fixed M x +, M^ 2 (and tan/3) one finds 
from eqs. (p.7| - |B.10|) that the remaining three neutralinos behave like M 2 . One then easily 
determines the (constant) behaviour of Mi from eq. ( |3.15| ) in the large M 2 limit, namely 
Mi ~ [M N2 (M 2 N2 - /x 2 - m|) - 4/x m| sin 2/3]/(M^ 2 - /x 2 ) 

5 Of course, more generally one could be interested in complex [i solution. The results of section 3.2 
and appendix B, not directly applicable in this case, would then need a generalization which is beyond 
the scope of the present paper. 
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In summary, the information from the plots in fig. 2 is rather interesting: apart 
from some small regions, the dependence of /i, M 2 (and even Mi to some extent) upon 
\M + — M + 1 is very simple for a wide range of the latter mass difference. In other words, 
once one chargino mass will be known with some accuracy (together with tan/?), fig. 2 
indicates that, in principle, the possible values of fi, M 2 , and Mi are strongly correlated. 
Moreover, although in fig 2 we have only varied M^+ the behaviour is quite generic and 
does not change much, qualitatively, when the other fixed inputs (Mn 2 and tan (3) are 
varied. 



Finally, as a last illustration of scenario Si, we plot in fig. ||] the resulting values of 
Mi and the three other neutralino masses, M Nl , M N3 , M N4 , determined from eqs. ( [3.161) - 



( |3.18|) . Note that the plots in fig. ^ are now functions of Mjv 2 , and with different M x +, 
M x + inputs than for fig. |[ The singularities of Mi for specific values of M N2 , correspond 
simply to the denominator of eq.( [3.15| ) vanishing (and also correspond to the exchange 
among neutralino masses, as illustrated on the plots). Apart from relatively small regions 
around the singularities, the almost linear dependence of Mi as a function of M Ni can be 
traced back to the fact that Mi behaves, for asymptotically large M N2 in eq.( |3.15| ), as 



Mi ~M Na + — [M Ni (M Ni - M 2 ) - / - M| + -Z- (/zM 2 - M^m2/3)] + 0(— r ) (4.1) 

Mn 2 M Ni N 2 



where is any of the neutralino masses other than Mjv 2 . Strictly speaking, ( [4.1| ) is 
valid only if the M^.'s become insensitive to the value of M^ 2 in this limit. That this is 
true can be easily seen from the fact that all the czj's in Eq.( p.7|) behave in this limit like 
~ M% 2 . Then the dependence on M N2 tends to cancel out in Eq.( |B.6| ), making the three 
neutralino masses insensitive to this parameter, as illustrated in fig. [3[ 



4.2 GUT scale values of the ino mass parameters 

Next, in fig. |] we illustrate plots very similar to the ones in fig. |2], but where in addition 
the parameters /x, M 2 , Mi have been evolved from a low scale, Qi ow ~ 200 GeV, to a GUT 
scale, Qgut ~ 2 x 10 17 GeV according to the renormalization group evolution (RGE) [|l"2]j . 
In order to illustrate as simply as possible some representative values after RG evolution, 
note that we have made a number of approximations for the evolution itself, that we 
feel are reasonable for the present purpose. Namely, the RG evolution is limited to the 
one-loop approximation, and we have assumed in addition a single, universal supersym- 
metric threshold, that is identified with the low energy scale, Qi ow ~ 200 GeV, where 
the evolution starts (so that the RG evolution is dictated by fully super symmetric beta 
function coefficients between Qi ow and Qgut)- There should be no difficulties in principle 
to incorporate in our framework more realistic supersymmetric threshold effects, provided 
however that the masses of other partners (Higgses and sfermions) are known at this stage 
in addition to ino masses. The latter refinements should however produce relatively small 
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corrections to the numerical values that are illustrated in fig. |j. 



It should not be very surprising that RGE scarcely change the shape of the various 
plots, with respect to those in fig. apart from a systematic shift and a slight change 
in the slopes. This is a direct consequence of the generic form of the one-loop evolution 
equations for the gaugino masses and the /x parameter, namely 



dm 
~dt 



where Qi refers to either gauge or Yukawa couplings and t is the evolution parameter. 
Taking into account the evolution of the couplings themselves, one then expects the mass 
parameters to behave like 

m(t)~m xn^W a 7£i(0) ai 

i 

where a, is some numerical power and mo the mass parameter at some initial scale. [The 
latter behaviour is exact, at the one- loop level, for the gauge couplings and also for the 
Yukawa couplings at least in the small tan (5 regime.] Thus the only effect of the running 
is a coupling dependent rescaling of the inversion results. 



4.3 Scenario S2 ' one chargino plus two neutralinos input 

Next we illustrate the (probably phenomenologically more relevant) scenario S2, namely 
with M x +, M N2 and M N , S as input. As expected, fig. || reflects the more algebraically 
involved inversion from the combined algorithm S" 2 (with unknown M x +), with the possible 
occurrence of several distinct solutions for (/x, Mi, M 2 ), as discussed previously. However, 
apart from the relatively untidy behaviour in some narrow zones, the domains of unique 
and twofold (or more) solutions are relatively well separated over a wide range of the 
Mn 2 values. Moreover, the behaviour of the plots for the particular input values in fig. |5] 
turns out to be quite generic. The shaded regions again corresponds to a zone where at 
least one of the output parameters (/1, Mi, M 2 ) becomes complex- valued. In fig. (H) we 
only show on purpose a range of values such that all masses are relatively light, while for 
larger M N2 the dependence of /1, M 2 , Mi upon the latter becomes simpler, almost linear, 
in accordance with the behaviour in the previous figures for scenario Si alone. There are 
also specific values of the input masses such that one of the solutions for /x, Mi, or M 2 is 
becoming very large, due to an explicit pole in the analytic inversion, as reflected e.g. in 
one of the /x and M 2 solutions, for M N2 ~ 80 GeV and M N2 ~ 123 GeV, respectively. 
Also, the dependence of scenario 5 2 upon tan/? variations is relatively mild. In contrast, 
varying M x + and/or M Nz input values can have more drastic effects, since in particular 
the occurrence of multi-fold solutions depends on those inputs. 

Finally, for completeness, we show in fig. ^ the values of the other physical masses 
M +, M^ and Mtv 4 , resulting from the same input choice as in fig. |5]. Note that, apart 
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from M Nl , the twofold valuedness of M+ 2 and M Ni for M N2 > 230 GeV turns out to be 
numerically negligible for this particular input. In a more complete and realistic analysis, 
it should be possible to examine e.g. whether some of the multi-fold solutions in fig. |5] 
could be excluded by looking at the consistency of the resulting other physical ino masses 
with data. 



5 Other soft breaking parameters inversion 

Ultimately, one would need to reconstruct the remaining part of the soft-breaking La- 
grangian, i.e. the sfermion and Higgs sectors. In contrast with the ino sector, the de- 
diagonalization of these sectors would not present much analytical difficulties (at least at 
the tree-level), provided one knows a sufficient number of physical masses. A full strategy 
(including loop effects) is, however, beyond the scope of the present paper. Our aim here 
is to simply sketch for completeness some straightforward inversion formulas and make 
some comments. 



5.1 Sfermion parameter inversion 

In the sfermion sector, de-diagonalization is straightforward since it involves 2x2 mass 
matrices if neglecting flavor non-diagonal terms. For instance, from eqs. ( |A.6|) , and as- 
suming the two physical stop masses m^, m f ~ 2 (with the convention < m^ 2 ) and the 
mixing angle d% to be given as input (plus tan/5), one immediately obtains 



A t = + (m~ - m~ ) (5.1) 

tan/5 V * 2 tlJ 2m t 



and 



Mq = m\ cos 2 9~ t + m| sin 2 9~ t - m\ - cos(2/3) (Am 2 w - m%) /6 (5.2) 

2 

M\ = m 2 x sin 2 9i + m 2 2 cos 2 9f — m 2 + - cos(2/3) (m 2 ^ — m|) 

Similar expressions are obtained for the sbottom and stau parameters. 
Evidently, in the case of a more constrained scenario, one would not need all sfermion 
physical masses and mixing angles to extract the relevant soft terms. In fact, even in 
the unconstrained MSSM case, Mq is a common parameter to stop and sbottom mass 
matrices , so that to reconstruct the complete soft terms of this sector one "only" needs 
five parameters, to be chosen among the four masses and two mixing angles. Conversely, 
notice that if knowing the six previous physical input, one can determine tanfi very simply 
from 

m^r cos 2(3 = cos 2 % + mj 2 sin 2 6i — mj- cos 2 Q~ h — mj| sin 2 Q~ h + m 2 — m 2 . (5.3) 
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5.2 Higgs sector inversion 



Going from physical masses to soft SUSY-breaking parameters in the Higgs sector also 
does not involve much difficulties, at least naively. In fact, electroweak symmetry breaking 
(EWSB) gives two (necessary but not sufficient) constraints, commonly used to fix the B 
and n parameter at the EWSB scale: 



m Hj - tan 2 3 m H ml 

" = Vi ^ ; (5 - 4) 

B\i = (rn 2 Hd + m 2 Hu + 2/i 2 ) sin 3 cos 3 



Now, one may assume for example that \x is extracted from the previous ino inversion, 
and use eqs. Q5.4p, together with the definition of the physical Higgs masses, eqs. ( |A.5| ), to 
determine the basic Higgs parameters. Typically, given the lightest CP-even higgs mass 
Mh, one obtains for the soft parameters at the tree- level: 



9 M 2 A — (u 2 + mi/2) (tan 2 3 — 1) 9 9 9 9 

m\ = — ^— — - ; mi = M 2 , — - 2a 2 5.5 

tf u 1 + tan 2 /? ' Hd A Hu P 

sin(2/3) 



5 = (m^ + m^ + W)^^ (5-6) 



and 



M fe 2 (m|-M 2 ) 
Ti| cos 2 2/3 - Ad 

and the other physical Higgs masses are determined from eqs. ( A. 5 ). Of course as is 



M 2 - hK z h > 7) 
Ma ~ m 2 z cos^20-M 2 



well known radiative corrections to the Higgs masses (in particular that of the lightest 
CP-even) are very important to take into account |13|. [T5| , for some values of the sfermion 
masses and mixing angles, so that it is much too unrealistic to restrict to such a tree-level 
Higgs parameter inversion. Including the leading one-loop corrections in expressions (|5.5| 
- 15.7| ) is in fact manageable. In this case the induced dependence on the full parameter set 
of the MSSM should be in principle taken into account. In particular, the parameters of 
the Higgs sector we are solving for enter now in a more complicated way equations ( |5.4| ) . 
However, a very good starting point would be to consider only the dominant top/stop- 
bottom/sbottom contributions in which case the corrections to Eqs.( [5.5| - |5.7| ) can still 
be written in an analytical form. Otherwise, one will still have to resort to a thoroughly 
numerical procedure. 



6 Discussion and outlook 

In this paper we have worked out a systematic inversion algorithm and strategy to ob- 
tain from relatively simple algebraic relations the basic MSSM Lagrangian, from some 
appropriate physical mass inputs. We mostly concentrated on the ino parameter inver- 
sion, which is the less straightforward algebraically, due to a complicated diagonalization 



17 



structure. The advantage over possibly more direct methods (like e.g. some systematic 
scanning of the Lagrangian parameter space), is that it gives a fast, reasonably tractable 
algorithm, and also clearly points out to the existence of non-trivial ambiguities in such 
a reconstruction, occurring in some definite range of the input mass values []. 
In addition, our results exhibit the relatively simple and generic behaviour of most of 
the Lagrangian parameters as functions of input masses, apart from some narrow regions 
where the behaviour is more involved. If developed further, performing more systematic 
simulations, this approach may thus potentially give useful insight into the precise connec- 
tion between the basic MSSM parameters and the experimentally measured ones. In the 
most optimistic case (that is, if knowing a sufficient number of physical input ino masses 
and with sufficient accuracy) it allows for a precise reconstruction of the unconstrained 
MSSM ino sector Lagrangian. In addition, some of the theoretically well-motivated ad- 
ditional assumptions, like for instance the universality of soft mass parameters, should 
be easily implementable in our basic algorithm, due to its relative flexibility. Typically, 
gaugino mass parameter universality at the GUT scale, leading to the approximate re- 
lation Mi(M z ) ^ 5/3 tan 2 6 , W M 2 (M Z ), can be very simply combined with our iterative 
algorithm 5*2, to give now /i, M 2 and M\ in terms of only one chargino and one neutralino 
mass. 



Finally we also mention the complementarity of our approach with the one discussed 
recently in ref. ||10|| . In particular, our simple analytical determination of Mi, eq. ( |3.15|) , 



may be readily implemented as well in their construction. Both procedures are however 
purely theoretical, in the sense that the influence of experimental errors on the input 
parameter measurements have not yet been taken into account. In any case, clearly an 
algebraic approach is a useful tool for a more systematic phenomenological study (taking 
into account e.g. the present bounds or the expected future accuracy on input parameters) 



to be presented elsewhere [TBI . 
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A MSSM mass formulas 

We collect here for completeness different useful expressions for the sfermion, ino and 
Higgs mass matrices or eigenmasses. 

-Gaugino sector: 
the chargino mass matrix reads from Eqs. (|2.4|, 127 



/ M 2 v/WsinM 
\ \J Im-w cos p jj, I 



3 Our main algorithm will be soon available as a computer fortran code, on request to the authors. 
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with the chargino mass eigenvalues: 



M\ x% = -[M 2 2 +/i 2 + 2m^±[(M 2 2 -/i 2 ) 2 +4m^cos 2 2/5+4m^(M 2 2 +/i 2 + 2M 2 /isin2/3)] 1/2 ] ; 

(A.2) 

and the diagonalization of the non-symmetric matrix ( A.l|) involves two mixing angles, 

tan20_ = 2VW »i Sin f + 5r S/ L (A.3) 



M|-// 2 -2m 2 v cos 2^ 

H cos /3 + M 2 sin f3 
' M%- ^ + 2m 2 w cos 2/? 



_ at /i cos /3 + M 2 sin f3 ,. A , 

tan 20+ = 2V2m WT7 £ — _ (A.4) 



-Higgs sector: 

the mass eigenvalues of the neutral CP-odd A, CP-even h, H and charged H + scalars 
read: 

M 2 A = m 2 Hu +m 2 Hd + 2i2 2 + A A ; (A. 5) 

M 2 H± = M 2 A + m 2 w + A H+ ; 

Kh = ^[M 2 +m| T [(Mi + m|) 2 -4M 2 m|cos 2 (2/?)] 1 / 2 ] + A^; 
where A a, Ah+, and A^h denote radiative correction contributions, whose complicated 



expressions in one-loop approximations are given e.g. in refs. 15, [L7 



-Sfermion sector: 
The mass matrices read 

M 2 = ( M 2 Q + m 2 + (lm 2 w -\m 2 z ) cos 2(3 m t ( A t - /// tan /?) \ 

1 \ rn t {A t — ///tan/3) m 2 + m 2 - |(m^, - m%) cos2/3 ) ^ ' ' 



M 2 = ( M^ + m 2 b - (\m 2 w + |m|)cos2/3 m b (A - /i tan/3) \ 

* y mb(A 6 — //tan/3) m 2 fl + m 2 + ^m 2 ^ - m|) cos 2/3 i ^ ' ' 

which after diagonalization give the stop and sbottom mass eigenvalues m^, m f ~ 2 ; mi , 
ra- h2 respectively, and their mixing angles % Q~ b . 

The f mass matrix has a similar structure as the b matrix, eq. ( A.7|) , with the replacements 
Mq — > Ml, mi — > mf K , m& — > m T , and appropriate D-terms. 



B De-diagonalization of M 

The following discussion can be carried out in terms of any pair of M^. (see section 3.2 
for notations). We choose for defmiteness M^^M^ and make the convenient relabeling 
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Mn 1 = X, Mtv 4 = Y and Z = XY. At this level we aim at the determination of some 
general consistency constraints and the distinction between input or output parameters 
is only for convenience. 



Using Eq. fl3.14D to substitute for Z in Eqs.( j3.1i| - j3.13| ), one obtains a set of three 
linear equations in, say, the variables X,Y,Mi. This system would give easily X,Y and 
Mi in terms of the other parameters in the form 



X 



A 



x 



Y 



A Y 



A 



Mi 



(b.i; 



A A A 

where it not for the fact that generically A = while Ax 7^ and Ay 7^ 0. Actually 
the same happens independently for the two subsystems Eqs.( |3.11] , |3.12j ) or Eqs. fl3.llT , 



when one tries to solve for X, Y, leading respectively to 



X 



Si 
5 



Y 



and 



X = — r ~ 




Y 



(B.2) 



(B.3) 



where 5 = and 5i, 62 are generically non vanishing. The necessary and sufficient 
conditions for the existence of solutions are thus Si = and 62 = 0. They can be cast in 
the following form: 



Pf j + (/1 2 + Mf - M X M 2 + (M x + M 2 )Sij - Sl)P tJ 
+//MK4M! + s 2 w M 2 ) sin 2(3 - [i 2 M x M 2 = 



(B.4) 



and 



( Ml + M 2 - S^Pl + (fi 2 (M 1 + M 2 ) + M|(c 2 M x + s 2 w M 2 - //sin2^))^- 
+/i(M|(c 2 u M 1 + s 2 w M 2 ) sin 2(3 - ^M 1 M 2 )S ij = 



(B.5) 



where Su = Mw. + Mw. and P,- 



M Ni M Nj , % = l,...,4,j = 1,..., 4 with i^j. Note 
that we used here the fact that equations ( |B.4| ) and ( p.5|) should hold for any pair of 
eigenvalues since we could have chosen as (A, Y) any set {M Nv M Nj ) with i ^ j. 

When Eq. (p.4j ) is satisfied, Eqs.( |3.1i|) and ( P-12j ) become equivalent to each other, 
so that the system made of ( |3.11| ), ( p,12| ) and ( 3.14]) is solvable in terms of X and Y. 
Similarly, Eq.( |B.5| ) does the same for the system ( 3.1 H ) ( p,13| ) and (|3.14 ). The bottom 
line here is that we have traded Eqs.( ^.ll| - |3.14p for the system ( |B.4| ), ( |B.5| ) which has, 
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however, some welcome features suitable for our purpose. An immediate consequence 
is that, [within the strategy defined in section 3 where M 2 and fi are determined from 
the physical chargino masses M x + , M x + and tan 0\ , one can determine uniquely Mi and 

any three neutralinos mass parameters, say M^ l , Mn 3 , Mjv 4 once the fourth, say M^ 2 is 
given. Indeed, substituting Mi from Eq.( [B.4| ) into Eq. ( |B.5| ) one obtains the following 
cubic equation in M^,i 7^ 2, 

a x M%, + a 2 M% + a 3 M Ni + a 4 = (B.6) 

where 



ai = M 3 N2 + M 2 {ii 2 -M 2 N2 )- M N . 2 + c 2 w M 2 z ) - c 2 wi iM 2 z sin 2(3 
a 2 = s 2 w M 2 z (M N2 - M 2 )(M N , 2 + fi sin 2/3) - M 2 a 1 

a 3 = s 2 w M 2 z (M% 2 + (M 2 - MivJ 2 // sin 2(3 + M 2 M Ni (M 2 - 2M N J) - {^ 2 + M 2 z )a x 
a 4 = fi[s 2 w M 2 z (M 2 (M 2 - M Na )(ji + M N2 sin2/3) + c 2 M| sin2/3(/isin2/3 + M^)) 
+(//M 2 -c 2 M! sin 2/3) ai] 



The three solutions of Eq. (|B.6|) are given by 



(B.7) 



M N3 = -L (_ a2+ (i+^_(lzi^ ) (B.9) 



3ax v ^45 2^2 
MAf4 -3^ ( - a2 + ^i^ 2^2^ ) (R10) 



with 



B= VC + V^ + C 2 (B.ll) 
A = -a§ + 3aia 3 (B.12) 



C = -2a2 + 9aia 2 a 3 - 27afa 4 (B.13) 

These solutions correspond necessarily to the three unknown neutralino mass param- 
eters Mnx-, Mn 3 , Mat 4 , due to the fact that, as we mentioned before, Eqs.( [B.4| , |B.5| ) should 



be satisfied for any pair (M^., Mjy.) with i 7^ j. Furthermore, injecting any of these solu- 
tions back into Eq. ( |B.4p determines the same and unique value for Mi given in Eq.( |3.15| ). 
Finally one can cast Eqs.( [B.^ - B.10| ) in an even simpler form once it is realized that 
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M Nl , M N3 , M Ni and M x are automatically real- valued when M 2 , /i (and M N2 ) are taken 
real. 

Indeed, when M N2 ,M 2 and fi are all real valued, i.e. ai,a 2 ,a 3 and a 4 real, the cubic 
equation ( |B.6| ) assumes at least one real solution. The real-valuedness of Mi follows then 
immediately from Eq.( 3.15| ). However, since the latter equation should be satisfied for 
any pair (M^., M^.) with i ^ j, then all the remaining neutralino mass eigenvalues must 
also be real. 



It thus follows that one can rewrite Eqs.(BT8 - B.1U) in an explicitly real valued form. 
To do this, we note first that the simultaneous reality of the Mjy.'s is equivalent to 



1 + ^4- 



A 
~B\'- 







(B.14) 



implying that A is necessarily a negative real- valued quantity. Plugging the above relation 
back into Eqs.Q - [BTT0|) one finds Eqs.(pl6| - CT) 
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Figure 2: /i, M 2 and M\ (with the "higgsino-like" convention < M 2 ) as functions of 
M x + for fixed M % + ( = 400 GeV); M Na ( = 50 GeV), and tan (3 ( = 2). The opposite 
"gaugino-like" situation, with > M 2 , is trivially obtained by the substitution M 2 — > \fj,\ 
and ji — > sign(/i)M 2 . 
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Figure 3: M\ and the three neutralino masses (Mn 1 , Mn 3 , Mtv 4 ) as functions of M^ 2 
fixed M + (= 80 GeV), M + (= 200 GeV) and tan/3 (= 2). 
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Figure 4: Same captions as for fig. 2 but with /x, M 2 and M 1 evolved up to GUT scale. 
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Figure 5: //, M 2 and M 1 as function of M° 2 = M N , 2 for fixed M Nz ( = -lOOGeV), M x + 
(= 80 GeV) and tan (3(= 2). 
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